Effective field theory approach to Casimir interactions on soft matter surfaces 
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We utilize an effective field theory approach to calculate Casimir interactions between objects 
bound to thermally fluctuating fluid surfaces or interfaces. This approach circumvents the compli- 
cated constraints imposed by such objects on the functional integration measure by reverting to a 
point particle representation. To capture the finite size effects, we perturb the Hamiltonian by AH 
that encapsulate the particles' response to external fields. AH. is systematically expanded in a series 
of terms, each of which scales homogeneously in the two power counting parameters: A = R/r , the 
ratio of the typical object size (R) to the typical distance between them (r), and 5 = ksT/k, where 
k is the modulus characterizing the surface energy. The coefficients of the terms in Art correspond 
to generalized polarizabilities and thus the formalism applies to rigid as well as deformable objects. 
Singularities induced by the point particle description can be dealt with using standard renormal- 
ization techniques. We first illustrate and verify our approach by re-deriving known pair forces 
between circular objects bound to films or membranes. To demonstrate its efficiency and versatility, 
we then derive a number of new results: The triplet interactions present in these systems, a higher 
order correction to the film interaction, and general scaling laws for the leading order interaction 
valid for objects of arbitrary shape and internal flexibility. 
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Objects which constrain a fluctuating field experience a 
Casimir interaction (l|. The underlying fluctuations can 
be either quantum mechanical or thermal in origin |2|. 
In this letter we will be interested in forces induced by 
thermal fluctuations between particles bound to surfaces 
characterized by surface tension (films) [3|-|6| or bending 
rigidity (membranes) |g-|9|. 

The non-trivial aspect of such calculations tends to 
arise from the constraints that the extended objects im- 
pose on the partition sum. This issue is usually dealt with 
by pinning the field to the surface of the objects through 
delta- functions in the integration measure [3[. A clear 
exposition of this method, applied to compact objects in 
fluid membranes and films, can be found in Ref. [6j]. For 
the electromagnetic Casimir effect this approach was re- 
cently improved and systematized within the framework 
of scattering theory, where the constraints of the objects 
enter the interaction energy through their scattering ma- 
trix coefficients [10(. 

In this paper we employ an effective field theory (EFT) 
formalism, orignally developed to study the gravity wave 
profile for inspiralling black holes [111, to streamline the 
boundary condition issue. This formalism has been uti- 
lized to derive not only new results in gravitational wave 
physics [12( but also to calculate the leading order fi- 
nite size correction to the Abraham-Dirac-Lorentz radi- 
ation reaction force law in classical electrodynamics [13Jj . 
Both of the these applications dealt with classical non- 
fluctuating fields, whereas here we will generalize the for- 
malism to allow for finite temperatures. 

To illustrate the effective field theory approach to soft 
matter surfaces most transparently, we will focus on con- 
straints imposed by mobile but rigid objects which pin 
field fluctuation modes along their circumference. Exten- 
sions towards more general types of constraints entails no 
change in formalism. 



We will assume that the bare surface Hamiltonian H 
is a quadratic functional of some field 4>(r): 



u[4>] = \k 



d 2 r <t>K<j) , 



(1) 



where k is a generalized modulus and the kernel K defines 
the physics of the problem. For films K = — V 2 while for 
membranes K = (V 2 ) 2 . In both cases 4>(r) is the surface 
height in Monge parametrization |l4j. 

The idea behind the EFT formalism is to treat the 
objects as point particles and to recapture their internal 
structure through additional terms Art = ^ Q C a O a in 
the Hamiltonian, where the scalars O a are polynomial 
in the field and its derivatives and a labels the parti- 
cles. The coefficients C a are chosen to reproduce the long 
wavelength physics, in analogy with block spin rcnormal- 
ization, and arc fixed by a matching calculation. In prin- 
ciple one must add all terms which are consistent with 
the underlying theory. On dimensional grounds the coef- 
ficients of these new terms will scale with powers of the 
particle size, R, such that the limit i? — >■ is well defined. 
Let us for the moment assume that the new terms obey 
a shift symmetry Arl(4>) = Art(4> + h), where ft. is a con- 
stant. This symmetry eliminates boundary conditions 
in which the object is pinned to a fixed height. Violat- 
ing this symmetry leads to long wavelength fluctuations 
which can lead to pathologies. Given this restriction, it is 
clear that we may truncate the sum in a derivative expan- 
sion, as each derivative will scale as A = R/r. Quadratic 
terms in <f> generate multipoles when the particles are 
subjected to external fields, hence we may interpret their 
coefficients as polarizabilities. Terms which are higher 
order in the fields will not have such a simple interpreta- 
tion. Such terms are suppressed by powers of 5 = k^T /k. 
To see this, we may simply rescale the field </> — > \/5 <j>. In 
this way the leading term in the partition function has a 



well defined <5 — s- limit and non-linearities are automat- 
ically suppressed by powers of 5. These ^-corrections will 
be treated in a forthcoming publication. Here we will be 
only concerned with the more canonical A-corrections. 
The leading order in A is unique and 0(A 2 ): 



AH^ 






Of* [dittrjdrfir*)] 



(2) 



The tensorial polarizability Cf allows for non-axisym- 
metric objects, but for the sake of simplicity we will re- 
strict to the symmetric case Cfj = C^Stj. The coefficient 
C® , which can be fixed by treating a single object in iso- 
lation, is the isotropic dipole polarizability. The reason 
for this nomenclature will become clear once we match 
for C®. Note that in the EFT formalism determining 
the proper Hamiltonian to reproduce the long distance 
physics is both conceptually and technically independent 
of finding the associated force between objects. We will 
therefore first discuss how to compute Casimir interac- 
tions between polarizable objects, as described by the 
Hamiltonian H c s = H + AH, and afterwards explain how 
the polarizabilities are fixed by a matching procedure. 

The canonical partition function of the system de- 
scribed by the effective Hamiltonian H c s is given by 



Z= f T><t> e-PK* 1 =Z (c 



-/3AM 



(3) 



where Zq is the partition function of the free Hamiltonian 
(AH = 0) and (• • • ) denotes the associated Gaussian 
average. Substituting AH from its definition in Eqn. (|2jl. 
one obtains the free energy of interaction, IA, through 



pu 



log(e-" A «) 




di<t>(r a )] 




(4) 



This cumulant expansion can be represented as a series 
of diagrams 15( , each depicting a (connected) 2n-point 



function of some derivatives of the field. 

The expansion (j4]) encodes n-body contributions to 
the free energy from the n th term onwards. For a given 
n-body force there are sublcading corrections stemming 
from higher multipoles as well as terms non-linear in the 
lower multipolc polarizabilities. 

Let us now illustrate this formalism by applying it to 
some cases of interest. We first consider circular particles 
on a surface with tension-dominated energy density, i.e. 
a film. The relevant differential operator is then —V 2 , 
with Green function G(r,r') = — j- log|r — r'\, and the 
modulus k will be replaced by the more familiar a for sur- 
face tension. At 0(A 4 ) we have the diagram in Fig. QJa), 
which represents the lowest order pair interaction. It 
Wick-contracts to 



pu (i) = 



E 

a . b 



4ff 2 



[didjG(r a ,r b )Y 



(5) 




FIG. 1: Relevant Feynman diagrams for the calculations in 
this letter. Solid lines correspond to propagators, as usual, 
while dashed lines represent world lines of particles. 



and evaluates to the pair potential 



w. 



(4) C a C b 

{ab} 4ir 2 o- 2 rj b 



(6) 



where we define r ao — r a — r b and denote by {. . .} an 
n-tuplet of distinct particles. Self energies corresponding 
to self-links in diagrams, such as the one in Fig. HJc), 
lead to divergent contributions to the interaction energy. 
However, for derivative interactions these divergences are 
all power like and can be absorbed into C®. These di- 
vergences carry no physical information, as there is no 
non-trivial renormalization group flow, and thus effec- 
tively we may set these diagrams to zero. 

Let us now consider contributions beyond 0(A 4 ) to the 
two body interaction, for which there are two possible 
sources: Higher order corrections to AH itself and higher 
order terms in the cumulant expansion. The next order 
term in AH, which is not sublcading in 5, is given by 



AH^^lj2 C ald t dMr a )Y 



(7) 



It corresponds to a quadrupole polarizability and is 
C(A 4 ). In principle there is also a term involving 
didi<f>(r a ), but it can be removed by a redefinition of 
4> [Ifili owing to the fact that di di<p = is the Eulcr- 
Lagrange equation for the problem. Inserting (J2J + © 
into the cumulant expansion Q generates a quadrupole- 
dipole interaction at C(A 6 ), given by 



/3< h} 



2(C a QC 6 D + C a D C b Q ) 



(8) 



At C(A 6 ) one might additionally expect a contribution 
involving three dipole interactions in the cumulant ex- 
pansion. However these vanish, since such terms neces- 
sarily involve self-energies, c.f. diagram [TJc). 

At C(A 8 ) we expect the quadrupole-quadrupole inter- 
action, arising from [d 4 G(r a , r b )] 2 : but also a non-linear 
dipole-dipole term, proportional to four dipole polariz- 
abilities and stemming from the 4 th cumulant, and dia- 
gram[U;e),with[a 2 G(r a ,r & )] 4 : 



(3U 



(8) 
{ab} 



7T 2 CT 2 327T 4 CT 4 



1 



(9) 



As long as the objects resist curvature, C^ ^= 0. However, 
if in addition to vertical translations the objects can also 
tilt, their ability to align with local gradients in </> implies 
C° = and the lowest order interaction is 0(A 8 ). While 
the precise value of C® might be hard to calculate for 
arbitrarily shaped rigid or elastic objects, Eqns. (|6I8I9[) 
nevertheless show that on films these will always interact 
with an asymptotic r~ s potential, thus generalizing a 
finding obtained in Ref. |5j for ellipsoidal objects. If the 
objects are identical, (C Q ) 2 > implies they attract. 

The lowest order 3-body interaction arises from three 
dipole insertions. Henceforth using the shorthand nota- 
tion G ab = G(r a ,ri,) and denoting partial derivatives by 
subscripts, the relevant interaction is 



P U {abc} 



^'a ^b c /~iab/~ibc /-tea 
2 U jfc fc? 



(10) 



which corresponds to diagramQJb) and yields a triplet in- 
teraction that scales as (r a ;,rb c r ca ) . However, owing to 
the symmetries of the tensor Gij — (8ij — 2fif.,)/(27rr 2 ), 
Tr[G n ] = if n is odd. Thus the leading dipole contribu- 
tion to any n-body interaction vanishes for n odd. The 
first non-vanishing 3-body interaction therefore arises at 
0(A 8 ), from terms which are quadratic in one of the 
dipole polarizabilities, diagram QJd): 



W\ll } 



167T 4 CT 4 
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4 4 
ab ac 



c» 



Cl 



'be 



4 4 
r cb r ca 



(11) 

This term drives triplets to attract, irrespective of their 
relative placement, as long as G D > 0. A possible dipole- 
quadrupole-dipole interaction also scales as 0(A 8 ), but 
again the symmetries of Gij and Gijki force it to vanish. 

So far we have calculated the forces in terms of a set 
of polarizability coefficients. These are fixed by calculat- 
ing an observable in the full (finite-sized particle) theory, 
expanding it in powers of A, and choosing the polariz- 
abilities such that the EFT reproduces the result to the 
appropriate order in A. The point to emphasize is that we 
are free to fix these coefficients by any means we choose. 
Thus it behooves us to choose as simple a setting as pos- 
sible: we will place a single particle in a simple stationary 
external field where we can easily calculate its response. 

To illustrate the procedure, let us match for C^ in the 
case of a rigid horizontal inclusion on a film. Consider 
placing the point particle in a background field such that 
the total field is given by 5(j)(r) + </>b g (r), where d(j)(r) 
is the induced field generated by the polarization of the 
inclusion. Terms linear in 54> in the Hamiltonian corre- 
spond to induced point sources 



p a (r) = -G a D ^ [6(r - r a )d^ hg (r)} 



(12) 



in the effective theory. The field emitted by an induced 
source is thus given by 

Mr) = -^ (a) G(r,r a )3f a Vb g (r a ) . (i 3 ) 



We may pick any background field we wish, but clearly 
it is simplest to choose the lowest multipolc field con- 
figuration necessary to generate a non-zero response. 
For a rigid horizontal inclusion this is a dipole field 
cj) h & ^ rcostp, i.e. one of constant slope. After solving 
the elementary boundary value problem (BVP) in the full 
theory and comparing to the effective theory result, one 
finds G D = 2irR 2 a. To match for G^, we need a back- 
ground with curvature, so we choose </> bg ~ r 2 cos(2iy9). 
Repeating the exercise we just performed, after appropri- 
ately generalizing (|T2"j) to account for the two derivative 
nature of 0, gives C Q = § i? 4 er. By Eqns. (|6l8l9l) this 
yields the pair interaction up to C(A 8 ): 



P U {ab} 



i? 1 



1 ab 



R 6 



1 <ib 



R H 



ab 



(14) 



The r^ 4 (dipole-dipole) and r^ 6 (dipole-quadrupole) 
interactions shown above agree with those derived in 
Ref. [4|. However, the r~ b 8 term consists of a lowest order 
quadrupole-quadrupole piece (prefactor "9" , also given in 
Ref. [J]) plus a non-linear dipole-dipole correction (pref- 
actor "1/2") that has not been previously calculated. Re- 
call that if the discs can also tilt, every term proportional 
to a dipole polarizability will vanish and thus only the 
interaction —9 k-oT{R/r a bf will survive. As discussed 
before, a r~ b -term will remain even for non-circular or 
bcndablc objects; only its prefactor will be different, ow- 
ing to the associated BVP being slightly different. 

Observe that the surface tension a (or generally the 
modulus k) cancels from the final result, because in ev- 
ery term of — f3U the number of Green functions, each 
accompanied by a factor a~ l , always matches the num- 
ber of polarizabilities, each proportional to a . This will 
no longer be the case once corrections beyond lowest or- 
der in 5 are included. 

Now consider the case of particles embedded in a sur- 
face with a bending-dominated energy density, i.e. a 
membrane. The kernel of the Hamiltonian in this case 
is K = (V 2 ) 2 , with the Green function G(r,r') = 
_L| r _ r '|2 l g| r — r'l, and the generic modulus k will be 
replaced by the bending rigidity k. We will assume that 
the particles can adjust to a constant slope background 
by tilting, so the first non- vanishing polarizability will be 
quadrupole in nature. Observe that this time, at 0(A 2 ), 
we need to write down two distinct terms: 



AH = 



sEl^w 



9^) 2 (r a ) + CQ'(a^^) 2 (r a ) 



(15) 
Terms involving the Laplacian cannot be removed via 
a field re-definition, because didi<fi = is no longer the 
Euler-Lagrange equation. Notice that due to the different 
kernel the quadrupole polarizabilities (7^ and C® scale 
~ R 2 , and not ~ R A as in the film case. 

The pair interaction follows easily from an expression 



analogous to Eqn. ([5]): 



-pu 



(4) 
{ab} 



c?c? 



2h : 



' (G ab 



ijkl 



, 2 i c a C b +C a C b ,^ ab 



2k? 



klj 



(16) 

The term ex C$ C^' {G$ kk f vanishes because Gf lkk = 
S(r a — r*b) by definition of the biharmonic Green func- 
tion. The other contractions are found to be (Gf^ kl ) 2 — 
4/(2wr 2 ab ) 2 and (G^ l f = 2/(2nr 2 ab ) 2 . 

The leading 3-body interaction in the membrane case 
stems from the third order cumulant. The interaction is 
derived from an expression analogous to Eqn. (fl"0| and is 
given by 



P U {abc} 



°a b °c [yijkl^klmn^r, 

\ ^ /^Q'/^rQ/^rQ I s~tab s-ibc si 
/ y °a °b °c [ U iikl < - r klmn ijr 



ca 
mnjj 



(17) 



where the sum includes three terms in which the prime 
accompanies a, b, or c. The full expression for this result 
is rather lengthy, but it simplifies greatly once we have 
the relationship between C® and C® . 

In the case of membranes we may perform a completely 
analogous matching procedure to the one described for 
the film. One finds the polarizabiliy coefficients C® = 
4ttR 2 k and C Q ' = -ttR 2 k, which through Eqn. dTHJ then 
yield the well-known pair potential [6J, |7| 



pu&> 



6 — ■ 

'ab 



(18) 



A similar force law could be derived for soft objects, after 
solving the relevant BVP to extract the polarizabilities. 
The potential would again decay ~ (R/r a b) , while its 
strength would depend on the relative elastic moduli of 
the objects and the bulk membrane, approaching 6 in the 
limit of rigid objects. 

With the polarizabilities fixed, we can simplify the 3- 
body interaction (|17p for the membrane, leaving 



^ {abc} 



-4 (c ab 



R 6 



Cbc 



-^ca) 2 2 2 
ab be ca 



(19) 



where c a b = cos[2(a a — ab)] and a a denotes the in- 
ner angle of the triangle at point a, etc. Unlike in the 
film-case, the sign of this triplet interaction depends on 
the relative orientation of the three particles. For in- 
stance, for an equilateral triangle of side-length s one 
finds W£? let (s) = +12 k B T(R/s) 6 , but for an isosceles 
triangle with angles (30°, 120°, 30°) and short side s one 
findsW t ( r 6 i ) plct ( S ) = -|fc B T(i?/ S ) 6 . 

To conclude, in this letter we have utlilized the effec- 
tive field theory approach introduced in [ll| to streamline 
calculations of Casimir forces on fluctuating two dimen- 
sional surfaces. We have reproduced well known results 
and derived several new ones, pertaining to non-linear 



corrections, 3-body terms, and general scaling laws for 
the leading order interaction between objects that are 
arbitrarily shaped and possibly flexible. However, this 
formalism extends well beyond the basic cases studied 
here and can be applied to many further situations in 
which existing techniques become rather unwieldy. For 
instance, the corrections beyond linear order in k^T /k 
follow in a straightforward manner and will be presented 
elsewhere. In addition, our results can be utilized to cal- 
culate the forces and torques between non-rigid objects 
or phase-segregated surface domains. This does not re- 
quire further field-theoretical sophistication but merely 
the calculation of their polarizability, which for compli- 
cated objects one might even decide to extract from ex- 
periment. 
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